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ABSTRACT 



Fourier analysis and power spectrum estimation of the cosmic microwave 
background anisotropy on an incompletely sampled sky developed by Gorski (1994) 
has been applied to the high-latitude portion of the 4-year COBE DMR 31.5, 53 and 
90 GHz sky maps. Likelihood analysis using newly constructed Galaxy cuts (extended 
beyond |6| = 20 deg to excise the known foreground emission) and simultaneously 
correcting for the faint high latitude galactic foreground emission is conducted on 
the DMR sky maps pixelized in both ecliptic and galactic coordinates. The Bayesian 
power spectrum estimation from the foreground corrected 4-year COBE DMR data 
renders n ~ 1.2 ± 0.3, and Qrms-PS ~ 15.3]![28 (projections of the two-parameter 
likelihood). These results are consistent with the Harrison-Zel'dovich n = 1 model of 
amplitude Qrms-PS ~ 18 /iK detected with significance exceeding 14cr {5Q/Q ^ 0.07). 
(A small power spectrum amplitude drop below the published 2-year results is 
predominantly due to the application of the new, extended Galaxy cuts.) 
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1. INTRODUCTION 

Following the seminal COBE DMR discovery of the anisotropy of the cosmic microwave 
background (CMB) after one year (1989/90) of observations (Smoot et al. 1992, Bennett et al. 
1992, Wright et al. 1992, Kogut et al. 1992) the instrument continued operating until late 1993. 
Intermediate, 2-year results of the mission were reported by Bennett et al. (1994) (also see e.g. 
Gorski et al. 1994, Wright et al. 1994). An overview of the 4-year results is provided by Bennett 
et al. (1996). In this paper we apply the CMB anisotropy power spectrum estimation technique 
introduced by Gorski (1994) to the definitive, 4-year 31.5, 53, and 90 GHz sky maps resulting 
from the COBE DMR mission. 

The COBE DMR-detected anisotropy of the CMB provides a unique opportunity to measure 
the spatial distribution of the inhomogeneities in the universe on the comoving scales ranging 

from a few hundred Mpc up to the present horizon size (inaccessible to any other astronomical 
observations) during the embryonic stage of their evolution, thereby avoiding the complications of 
(cosmologically) recent evolution of most astrophysical systems. There are, however, important 
experimental/observational limitations on the extent to which cosmologically interesting quantities 
can be extracted from the DMR data. These are imposed by our understanding of systematic 
effects, instrTimcntal noise properties, and non-cosmological foreground signals peculiar to our 
location in the Galaxy and in the universe, and were investigated as follows. Kogut et al. (1996c) 
have addressed systematic effects in considerable detail and found them to be unimportant relative 
to the accuracy to which cosmological parameters can be determined from the DMR 4-year data. 
Gorski et al. (1996) provide a comprehensive description of the noise properties of the DMR 
sky maps and demonstrate the sufficiency of the DMR noise models (see also Lineweaver et al. 
1994). Emission from the galactic plane cannot yet be modelled to adequate precision and must 
be excised from the data (see §2 and Banday et al. 1996a). Foreground galactic emission at 
high latitude (Bennett et al. 1992, Kogut et al. 1996a, b) can be partially accoTintcd for by using 
Galaxy-dominated spatial templates at frequencies where the CMB does not dominate (see §3, and 
Kogut et al. 1996a,b). The potential contamination of the CMB anisotropy by the astrophysical 
foregrounds outside of our Galaxy is demonstrated to be negligible in Banday et al. (1996b). 

In this work we focus on a general characterization of the 4-year DMR sky maps in terms of a 
power spectrum of the angular distribution of the CMB temperature fluctuations and attempt to 

estimate its sensitivity to a plausible foreground emission contamination. We follow and further 
develop the method introduced by Gorski (1994, hereafter G94a) and applied to the 2-year DMR 
data in Gorski et al. (1994, hereafter G94b). As a useful power spectrum parametrization, we 
utilize a power-law family specified by Qrms-PS and n — the amplitude and shape parameters 
(Bond & Efstathiou 1987, Fabbri, Lucchin, Sz Matarrese 1987), but we abbreviate the symbol 
Qrms-PS and use Q instead. The results of the analysis involving the CMB spectra of specific 
large scale structure models will be presented elsewhere. Hereafter, bold letters denote matrices 
(upper case) and vectors (lower case), and p is a pixel label. 
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2. DATA 

Our published analysis of the 2-year DMR data (G94b) involved a simple |6| > 20° Galaxy 
cut, which retained the relatively bright foreground emission from the Scorpius-Ophiucus and 
Taurus-Orion regions (de Vaucouleurs 1955). In this work we employ new, extended galactic cuts 
(see Banday et al. 1996a) which further excise these regions from the data used for the inference 
of the cosmological CMB anisotropy. We consider both extended |6| > 20° (denoted '20+') and 
extended |6| > 30° ('30-1-') Galaxy cuts. We analyze the DMR sky maps in both ecliptic ('E') 
and galactic ('G') pixelization. To facilitate a comparison of the 4-year results with the published 
analyses of the 2-ycar data we re-apply the simple \h\ > 20° Galaxy cut. From a total of 6144 data 
pixels on the entire sky, the following numbers of pixels remain for analysis after applying the 
different galactic cuts: 4038 for E20, 4016 for G20, 3890 for E20+, 3881 for G20+, 3039 for E30+, 
and 3056 for G30+. 

We form weighted average maps at each frequency by combining the A and B channels, 

A31 = 0.611 A3L4 + 0.389 A31B, 

A53 = 0.579 A53A + 0.421 A53B , (1) 
Ago = 0.382 A90A + 0.618 AgoB, 

minimizing the cut- sky- averaged noise variance per pixel in the combined maps. The weighted 
31.5, 53, and 90 GHz 4-ycar cut-sky DMR maps are analyzed jointly (equivalent to coaddition) 
with and without the foreground correction, or linearly combined (after foreground correction, 
indicated by tilde) with the weights designed to remove any leftover free-free emission: 

AcMB = -0.302 A31 + 0.633 A53 + 0.669 Aqq. (2) 

The weights given to both 31.5 and 90 GHz maps in this linear combination technique are 
considerably larger than those in the coadded technique (~ 0.08 and ~ 0.26, respectively). As a 
result, the noise is larger in the linear combination map Aqmb, slightly exceeding the rms noise 
in the 90 GHz map alone. 

Fourier analysis of the cut-sky maps is performed (see G94a for details) in the 961-dimensional 
linear space spanned by the orthonormal functions ^, which are linear combinations of spherical 
harmonics with i < 30. A sky map at frequency v, Ai,{p) (where p is the pixel label) is Fourier 
decomposed as follows 

^''"^6144: ^ ^i^i.P)Mp) and A^{p)\p^^^^t .^y} = ■ ilf{p). (3) 

p&{cut sky} 

Fourier coefficients Ci are linear in pixel temperatures, and account for the cut-sky mode coupling 
explicitly. Figure 1 (Color plate) displays an example of the ^ basis function constructed for 
various versions of the Galaxy cut used in this work. 
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The dominant component (in rms per pixel) of the raw, dipole-subtracted, cut-sky DMR maps 
is the instrument generated noise. A correct statistical modeling of this contaminant is essential 
for any method attempting to quantify the cosmological CMB anisotropy. A detailed analysis of 
noise properties of the entire DMR data set is discussed in Gorski et al. (1996). Here we employ 
the noise modeling in Fourier space as described in G94a,b. The noise content of the 4-year DMR 
data is characterized by the rms noise amplitudes per Fourier mode of approximately 10.5, 3.4, 
and 5.5 iJ,K at 31.5, 53, and 90 GHz, respectively. The coaddition of all available measurements at 
three frequencies renders a map with the rms noise amplitude per Fourier mode of ~ 2.7 fiK. 



LIKELIHOOD ANALYSIS OF FOREGROUND CORRECTED DMR SKY 

MAPS 



Following G94a,b, we conduct Bayesian inference from the COBE DMR data of the power 
spectrum of primordial CMB anisotropy under the usual assumption that the measurement (here 
the joint vector of coefficients of Fourier decomposition of the 31.5, 53, and 90 GHz sky maps, 
= (cg'j, Cg^g, Cq^q)) Is a sum of a signal — primordial CMB anisotropy, and noise — a random, 
receiver generated contamination, and that both these components are Gaussian stochastic 
variables. Hence, the likelihood function is exactly represented as a multivariate Gaussian in 
variable Cy. The likelihood density is fully specified by the Fourier space correlation matrices 
of signal, Cs, a function of the assumed (and tested) power spectrum of the CMB anisotropy, 
and noise (at a given frequency u), Cjv^,, which are are computed as described in G94a,b. No 
approximations or Monte Carlo simulations are required for evaluation of the likelihood function, 
and this power spectrum parameter estimation technique, bi-linear in the measurement values of 
temperature perturbations, is statistically unbiased. 

This formulation of the likelihood problem allows one naturally to address a problem of 
contamination of the CMB anisotropy data by the diffuse foreground emission. Assume (generally) 
that the anisotropy, Cj/, is measured at M frequencies v, and the foreground emission can be 
measured or modeled as a certain number, N, of non-stochastic spatial templates Fk{p) (whose 
Fourier transforms are ffe, k = 1,. . . ,N). The template corrected joint data vector, Ctu, and its 
covariance matrix, Csjvtti, are 



I 



Ctti 



s-^k=h ft e 



I. r. 

a 



Csivtd = (ctti- Cy ) = 



where the matrix C 



C5 + C^f^ 



C5 



(4) 



SNy- = + specifics the Gaussian probability distribution of the 

receiver noise contaminated theoretical CMB anisotropy signal at a given frequency Vi. With these 
definitions the likelihood function takes the usual form 



P(ctij) oc exp(-cj • Cg^y • ctij/ 2) / v^detQ^. 



(5) 
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This is a function of all parameters used to describe the CMB power spectrum and M x N 
coefficients used in the foreground correction. For a given CMB power spectrum (i.e. a fixed 
C5) the likelihood can be maximized with respect to parameters usual fitting problem 

(naturally, the dimensionality of the data vector is presumed to be large compared to both M and 
N). Prom the maximum likelihood equations, d InP/dal = 0, one derives the following set of 
linear equations (for a given u and k): 

E C • (csk^) ■ f. = E C • (csk^) ■ e." , (6) 

j,u' v" 

which render the maximum likelihood solutions S^. The parameter covariance matrix (dimension 
(M X Nf) is 

In both equations (c^^^^ ^ is a square portion of the inverse covariance matrix C^^^ acting 
between frequencies and v' . 

We apply this formalism to study the sensitivity of parametrization of the cosmological 
anisotropy derived from the DMR data to the possible contamination by high latitude galactic 
foreground emission. We assume that the faint foreground glow of synchrotron, free-free, and dust 
emission outside the galactic plane can be traced by the available templates from observations at 
frequencies where the CMB docs not dominate. Two plausible selections (Kogut et al. 1996a, b) 
are the DIRBE 140 jiva data (Reach et al. 1995), and the 408 MHz full sky radio map (Haslam 
et al. 1980). (We also used the 110, and 240 /xm DIRBE data and found the coupling of these 
to the DMR frequencies indistiguishable from that derived at 140 yum.) We conduct the analysis 
as outlined above with M = 3 DMR frequencies and N = 2 templates. Together with the CMB 
parameters Q and n this constitutes an 8-dimensional likelihood problem. (A further extension of 
the method, described in Banday et al. 1996b, addresses the issue of cross-correlation of the DMR 
data with existing catalogues of extragalactic objects; there we have employed N = 6 templates, 
and solved the 20-dimensional likelihood problem.) For each Galaxy cut selection we evaluate the 
Q-n likelihood both without and with the foreground correction applied to the data. At each Q-n 
grid point we derive the maximum likelihood solutions a^(Q, n) (and their covariance matrix), and 
use these values to evaluate the foreground corrected likelihood. Maximization of the corrected 
likelihood renders simultaneously the most likely solution of both the cosmological (Q-n), and 
foreground (a'l) fitting problem. The likelihood surface in Q and n is sharply peaked (as shown in 
G94b for 2-year data), and despite it being non-Gaussian in functional form, the average values 
of the foreground coupling constants {{al{Q,n)), weighted by the corrected likelihood) turn out 
practically identical to the maximum likelihood solutions. The likelihood averaged values of the 
six coupling constants between the three DMR maps and the 140 fim and 408 MHz data are shown 
in Fig. 2. The physical interpretation of these results in terms of the inferred contribution of 
galactic dust, free-free, and synchrotron emission in the high-latitude DMR sky maps is discussed 
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by Kogut et al. (1996a,b). Here, we proceed to discuss the primary target of the analysis of the 
foreground corrected DMR sky maps — the cosmic CMB anisotropy. 



4. RESULTS 

Fig. 3 (Color plate) shows the coadded, foreground corrected DMR sky maps in both the 
ecliptic and galactic coordinate pixelization. For visualization the maps have been "cleaned" by 
removing the £ > 40 noise (the 7° FWHM DMR beam picks up sky signals that are entirely 
contained within £ < 40 in Fourier space). Ecliptic and galactic maps appear very consistent at 
large angular scales, but near the pixel scale one can notice the localised effects of different noise 
binning depending on the reference frame. The choice of pixelization re-bins roughly half of the 
instrument noise while leaving the CMB anisotropy unchanged; consequently, slight variations in 
the likelihood results are expected. 

Fig. 4 displays the cut sky DMR power spectra before and after the foreground correction 
evaluated for four combinations of pixelization and Galaxy cut. (These spectra are shown here for 
illustration purposes only; the likelihood method uses the Fourier coefficients q, not any quadratic 
function thereof.) The apparently steep power spectrum of the foreground Galaxy templates, 
~ at higher was discussed by Kogut et al. (1996a). An important result of the simultaneous 
CMB background {Q, n) and foreground (a^) fitting is that the Galaxy correction affects only the 
very low order multipoles of the derived cosmic CMB anisotropy. Since the quadrupole moment 
of the Galaxy is partially counter-aligned to the CMB quadrupole (Kogut et al. 1996b), correcting 
the DMR maps for Galactic emission increases the estimated CMB quadrupole. 

Table 1 provides a comprehensive summary of the maximum likelihood solutions for the CMB 
anisotropy power spectrum model parameters. The first entries in Table 1 refer to the high-latitude 
coadded DMR maps from which no model of Galactic emission has been removed and include, for 
comparison with the published 2-year analyses, the results of re-application of the \h\ = 20 deg 
Galaxy cut. One should note the following: (1) compared to the G20 2-year results described in 
G94b, the 4-year DMR data render a slightly smaller, but statistically consistent, power sepctrum 
normalization (roughly a ~ 0.5 jiK drop of Q|n=i); (2) solutions for n are fairly stable with respect 
to variation of the Galaxy cut, but the derived anisotropy amplitude decreases by more than 1 ^uK 
upon removal of the bright Scorpius-Ophiucus and Taurus-Orion foreground regions (i.e. upon 
extension of Galaxy cut from '20' to '20-|-'), (3) an extra, but smaller, decrese of the amplitude 
results from a further extension of Galaxy cut (to '30+'); we consider this as suggestive of having 
reached a sufficiently deep Galaxy plane excision with the '20-I-' cut. Certainly the higher Galactic 
cut excludes real CMB features correlated over large areas of the sky in the range 20° < |6| < 30°, 
which affects the Q — n fits. 

Comparing results for which the quadrupole is included {imin = 2) to those from which it 
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has been excluded {£min = 3), we find again a systematic shift to higher normahzation and lower 
index n in the quadrupole-excluded results. The shift is well within the uncertainties from noise 
and cosmic variance, but reflects at least in part the counter-alignment of the CMB and Galactic 
quadrupoles: if the quadrupole is to be retained in the inference of cosmological parameters, 
Galactic emission should be accounted for in the data. 

The next set of entries in Table 1 refers to the DMR maps corrected for Galactic emission 
using the 408 MHz survey to trace synchrotron emission and the DIRBE 140 fim map to trace dust 
emission. Free-free emission is removed using a linear combination of the dust- and synchrotron- 
corrected DMR maps. This map makes no assumptions as to the distribution of free-free emission, 
but is relatively noisy compared to the coadded maps, which is reflected in lower significance of all 
the derived parameters. The final set of entries refers to the maximally sensitive coadded maps, 
corrected for Galactic emission using the 408 MHz survey to trace synchrotron emission and the 
DIRBE 140 //m map to trace both dust and free-free emission. The fitted values for Q and n 
are in good agreement for both foreground correction techniques. One will note that the high 
latitude foreground correction further reduces the fitted power spectrum amplitude by ~ 0.5 fxK. 
This cannot be considered as significant particularly in the light of a similar order of discrepancy 
between the galactic and ecliptic map maximum likelihood solutions, and especially in comparison 
with an overall statistical significance of all the determined parameters. 



5. SUMMARY 

A statistically unbiased technique of power spectrum determination has been used to infer 
the CMB normalization Qrms~PS s^nd power-law index n from the COBE DMR 4-year sky maps. 
The sky maps are decomposed in the basis of orthonormal functions on the cut sky into a set 
of coefficients linear in the pixel temperatures. These coefficients are then used in a likelihood 
analysis to infer the parameters of the CMB anisotropy models. We account for high-latitude 
Galactic emission by using the 408 MHz sky survey to trace the angular distribTition of synchrotron 
emission and the DIRBE 140 ^va. map to trace dust emission. Free-free emission is removed using 
either a linear combination of the DMR maps or using the DIRBE map to remove free-free as well 
as dust emission. Since we find some minor changes in the fitted parameters with the Galaxy cut 
or map pixelization and method of high latitude foreground modeling, a judgement call is required 
to quote a unique COBE DMR determination of the primordial power spectrum parameters. The 
4-year data are consistent with, qucidrnpole normalizcition Qrms— PS ~ IS.StIg ijK and power-law 
index n ~ 1.2 it 0.3. With n fixed at unity we derive the normalization for an exact scale-invariant 
Harrison-Zel'dovich power spectrum Qrms-Ps\n=i ^ 18 /uK with high ( ^ 14a) significance 
6Q/Q ^ 0.07. Implications of the COBE DMR 4-year data for the specific models of large scale 
structure formation will be discussed elsewhere. 
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FIGURE CAPTIONS 



Fig. 1 (Color Plate). An example of the full-sky spherical harmonic and the related cut-sky 
orthonormal functions (Gorski 1994) — note the coordinate frame and Galaxy cut dependence — 
used in Fourier analysis of the DMR and foreground emission template data. 



Fig. 2. Coefficients of coupling between the three DMR channels at 31.5, 53, and 90 GHz and 
the Galactic templates: the 140/im DIRBE data — filled symbols, and the 408 MHz radio data 
— open symbols. The coefficients shown in the plot (computed for the choice of pixelization and 
Galaxy cut as indicated) are the maximum likelihood values averaged over the {Q — n) likelihood 
of the simultaneous background and foreground fit to the 4-year DMR data. The one a confidence 
intervals are defined in the usual manner as the square roots of the diagonal on the covariance 
matrix. 



Fig. 3 (Color Plate). COBEBMR 4-year, 31.5, 53, and 90 GHz sky maps coadded after foreground 
correction according to the 140/im- and 408 MHz-template fitting. Galactic and ecliptic (rotated 
to galactic coordinates to facilitate the comparison) pixelizations are shown. Noise reducing 
Fourier removal of the i > 40 content of the raw data has been applied to the full sky maps in 
order to aid visualization. Apparent small differences in pixel temperatures are induced by the 
reference frame dependent time ordered data noise binning in the map making procedure. 



Fig. 4. Power spectra of the coadded DMR data and the templates used to assess the 
galactic foreground emission. The four panels exhibit the cut-sky rms power spectra, 

C£ = (^IZp^i c?^ / {2i + 1) , evaluated for both ecliptic and galactic data with two Galaxy 
cuts as indicated. (Small letter q is used in the definition of the cut-sky spectrum to avoid 
confusion with commonly used Q symbols, which denote the theoreticall, full-sky anisotropy 
spectrum.) Open circles — coadded raw DMR data, filled circles — template corrected DMR 
data, filled squares — foreground emission contribution to the DMR 53 GHz data characterized 
by spatial structure of the 140^m DIRBE map, arrows — 95% upper limit to possible foreground 
emission in the DMR 53 GHz channel characterized by spatial structure of the 408 MHz radio 
data, horizontal lines — rms noise per Fourier mode in the coadded 4-year DMR maps. 



Table 1: Derived Parameters from the Likelihood Analysis of COBE DMJi Four- Year Sky Maps' 
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Linear Combination Map 
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Coadded Map: galaxy corrected 
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-2.83 -5.08 

,7 +1.39 +2.97 

-1.34 -2.59 



We derive the tabulated values for Q and n by projecting the 2-dimensional, Q — n likelihood 
function onto the Q ov n axis (rather than marginalising) and quote the maximal likelihood 
value and the 68% and 95% credible intervals of the resulting 1-dimensional likelihood 
densities. We also evaluate the likelihood along the slice n = 1 to derive the CMB 
normalization for an exact Harrison-Zel'dovich power spectrum. 
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